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We consider a few-boson system confined to one dimension with a single distinguishable particle 
of lesser mass. All particle interactions are modeled with (5-functions, but due to the mass imbalance 
the problem is nonintegrable. Universal few-body binding energies, atom-dimer and atom-trimer 
scattering lengths are all calculated in terms of two parameters, namely the mass ratio: m^/mn, 
and ratio <?hh/</hl of the 5-function couplings. We specifically identify the values of these ratios for 
which the atom-dimer or atom-trimer scattering lengths vanish or diverge. We identify regions in 
this parameter space in which various few-body inelastic process become energetically allowed. In 
the Tonks-Girardeau limit (jhh —> oo), our results are relevant to experiments involving trapped 
fermions with an impurity atom. 


I. INTRODUCTION 

Strongly interacting one dimensional (ID) quantum 
systems have been of fundamental interest for many 
years mm- Not long ago, a number of experiments 0- 
[5] involving quantum gases tightly confined to two- 
dimensional (2D) optical lattices realized the ID Lieb- 
Liniger-McGuire model EHS]- Interpretation of these ex¬ 
periments in terms of the one-dimensional model param¬ 
eters has been facilitated by a calculation of the effec¬ 
tive ID coupling constant (i?id) and dimer energy (E 2 ) 
in terms of the three-dimensional (3D) s-wave scatter¬ 
ing length ( 030 ) and the transverse confinement length 
(aj_ = (^w_ l) -1/2 ) I3EEQ] (in units with h = 1 through¬ 
out): 

_ 2a 3D * /1 n _ aj - ri\ 

510 Mai(l-|CQ)|ff)’ a 3D (1) 

where £ (z,q) is the generalized Riemann zeta func¬ 
tion HU, p. is the two-body reduced mass, and UIb = 
1/2 — E 2 /(2uj±) is the dimensionless dimer binding en¬ 
ergy. The result predicts the so-called “confinement 
induced resonance” (CIR) when |£(l/ 2 )|a 3 D = aj_, 
where |(/id| —> 00 and the ID scattering length aiD = 
— (/t^id ) -1 passes through zero. Remarkably, g id can be 
experimentally tuned by varying either a 3 D (via a mag¬ 
netic Feslibach resonance), or aj_. 

Moreover, equation 0 predicts that a ID dimer al¬ 
ways exists below the asymptotic threshold energy lo± 
regardless of the sign of 030 - Such confinement-induced 
dimers have been seen in Monte Carlo simulations 
and observed experimentally by RF spectroscopy m- 
When the axial extent of the dimer is large compared to 
the transverse confinement length, the few-body physics 
is expected to follow from a purely ID calculation with 
(5-function couplings given by Eq. 0 pu m na nsj. 
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This limit is achieved when aj_/a 3 D —> — 00 , such that 
UIb 1. 

More recent quasi-ID experiments have begun to probe 
the interaction of a degenerate quantum gas with a distin¬ 
guishable impurity mmm- Motivated by the afore¬ 
mentioned evidence that—in the appropriate limit— 
purely ID calculations of few-body processes can pro¬ 
vide physically meaningful insights into quasi-ID exper¬ 
iments, we consider here the N = 3 and N = 4 instances 
of the following Hamiltonian: 

H _ _(L v — _ A_ — 

2toh ^ dxj 2?7?l dx\ 

(2) 

TV—1 TV-1 v ' 

+ X 9tmS(xi - Xj) + ^2 gm,S(xi - x N ), 

i<j i=1 

where particles 1 through (N— 1) are identical bosons (H) 
of mass mu, and particle N is an impurity (L) of mass 
TOl- The equal-mass (tol = mn) instance of Eq. ([ 2 ]) 
has been realized with 40 K atoms [14] and more recently 
6 Li [T7|, with related theory work found in Refs. 

Note that the Bose-Fermi mapping [251 allows one to con¬ 
sider the (?hh —> 00 instance of Eq. pi for a description 
of a noninteracting background gas of fermions. The un¬ 
equal mass case (mn > m l) of Eq. 0 has been realized 
with a 41 K impurity in a background gas of 87 Rb [15] , 
We assume that the impurity is of lesser mass than the 
background atoms, and employ the Born-Oppenheimer 
method to obtain adiabatic potentials for the heavy- 
particle motion. In the thermodynamic limit, stabil¬ 
ity of the background gas requires jhh > 0. At the 
few-body level, our calculation permits us to consider 
<7 hh < 0, but with the caveat that the lowest scatter¬ 
ing threshold involves a bound impurity. The eigenstates 
of Eq. pi are completely specified by the coupling ratio 
A = <?hh/9hl and the mass ratio /? = TOl/toh- The case 
N = 3 of Eq. 0 has been studied in considerable detail 
by Kartavtsev et al [3S|, who treat all possible values of A 
and P using the adiabatic hyperspherical representation. 
Mehta [30 a used the Born-Oppenheimer method to ex¬ 
tend those results to the case N = 4, but only considered 
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the cases A —► 0 and |A| —► oo with /? < 1. This work 
is an extension of those calculations to finite values of A. 
Some of the approximate N = 3 calculations presented 
here can be directly compared to the benchmark calcula¬ 
tions of |29j . but we are unaware of similar calculations 
for N = 4. 

We calculate the atom-dimer scattering length oad, 
the atom-trimer scattering length oat and trirner and 
tetramer bound-state energies as a function of /? and A. 
The critical values of /3 and A at which these scattering 
lengths diverge are identified and marked by the appear¬ 
ance of a new bound state. We find that few-body bound 
states tend towards deeper binding as one increases A 
(making it less negative) or increases /3 _1 . This behav¬ 
ior allows us to identify the critical values of A and /? 
where particular inelastic processes become energetically 
allowed. 


II. BORN-OPPENHEIMER SOLUTION 


Here, we briefly sketch the Born-Oppenheimer calcu¬ 
lation. Many elements of the derivation are unchanged 
from ED], so we present an abridged derivation sufficient 
to highlight changes due to choosing finite A = <?hh/3hl- 
We scale the Hamiltonian by the HL binding energy, 
B 2 = /xhlShl/ 2 = l/(2/i HL a 2 ). and use the HL scat¬ 
tering length a = — 1/(/xhlShl) > 0 as the fundamental 
length unit. The calculation of few-body observables in a 
homogeneous system begins the removal of the center-of- 
mass motion by a transformation to mass-scaled Jacobi 
coordinates: {x,y,z}:. 



- x 3 



94b V TO i + m 2 

Mi23,4 ( m 1 xi + m 2 x 2 + m 3 x 3 _ A 
94b V m i+m. 2 + m 3 ) ' 


(3) 


Here, p 4 b = ( 912912 , 39123 , 4) 1 / 3 is the four-body reduced 
mass, and p a ,b is a two-body reduced mass between clus¬ 
ters a and b. The coordinates for the three-particle prob¬ 
lem are the same, but with the ^-coordinate omitted and 
94 b replaced by p 3b = (pi 2 yi 2 , 3 ) 1/2 - 


A. Three-Body Problem 

The Born-Oppenheimer factorization of the three- 
body wavefunction ty(x, y) = &(x',y)ip(x) requires that 
<l>(x;y) be an eigenstate of the fixed-x Hamiltonian: 

H ad =-^ L3 ^ x 2+93[S(y-x 0 ) + S{y + x 0 )\, (4) 

with x-dependent eigenvalue u(x). We have defined 
the unitless three-body reduced mass p 3 = (1 + 


/3)/[2 \//3(2+ /?)], coupling g 3 = -2\/2[/?/(2 + /3)] 1/4 , and 
scaled heavy-particle separation xq = X\ff}/(2 + /?). The 
solution to Eq. @ yields the following transcendental 
equation for the lowest u(x): 

—--b 1 = — exp(— 2kxq), (5) 

9393 

where n = — 2pm J u(x). The heavy-particle wavefunc¬ 
tion ip(x) is now governed by the effective Hamiltonian: 

= 2p 3 h 2 + 93X6 ^ + < X) + (6) 

Here, Q(x) = (<&' |$') is the diagonal nonadiabatic cor¬ 
rection, where the primes denote derivatives with respect 
to the slow coordinate x , and the integration is carried 
out over the fast coordinate y only. Ignoring Q(x) leads 
to a lower bound on the three-body ground-state energy, 
while including Q(x) leads to an upper bound that is typi¬ 
cally more accurate [5TH33] . We refer to the former calcu¬ 
lation as the “extreme adiabatic approximation” (EAA) 
and the latter as the “uncoupled adiabatic approxima¬ 
tion” (UAA). The cases A = 0 and |A| —> 00 were treated 
hr [3D]. Here we consider arbitrary A, and replace the 6- 
function in Eq. Q with the following boundary condition 
on the solution ip(x) for x > 0: 

A 1 + & m 
\ ipdx ) x ^ 0+ v ^0[/3(2 + /3)] 1 /4 


B. Four-Body Problem 

For the four-particle problem, it is convenient to work 
in “hyper-cylindrical” coordinates, trading { x , y, z} —> 
{p,(j),z} by the usual transformation: p 2 — x 2 + y 2 , and 
tan 4> = y/x. We make the Born-Oppenheimer factoriza¬ 
tion: ^(p, </>, z) = $(p, (/); z)ip(p , (/>), and integrate out the 
light-particle degree of freedom by demanding <E>(p, </>', z) 
be an eigenstate of the following fixed-{p, </>} Hamiltonian 
with eigenvalue U ( p , <j>): 

H ^ = J^^ +9A ^ S{Z ~ Zi) - (8) 

r i =1 

We have again introduced unitless parameters: P 4 = 
(P + l)/[2/3 2 / 3 (3 + /3) 1 / 3 ], and y 4 = -2^3 \p/(3 + /3)] 1/3 . 
T he heavy-li ght coalescence points occur at z t = 
y/20/(3 + /?) psm(cj>— fa) with <j> 1 = —in/3, 4 >2 = 0, 
and <p 3 = —2tt/3. The adiabatic equation is simply a 
triple-<5 problem that leads to the following transcenden¬ 
tal equation for the potential energy surface: 

(9 d~ 2/4.) ^ 2 kz 3 _|_ (g ~ 2^) c 2kzi _ e 2KZ 2 _|_ e 2K(zi — z 2 + z 3 ) 

g 2 9 + 2k 

(9) 

where g = 2p 4 y 4 , k 2 = — 2p 4 /7(p, (j>) > 0. 
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The heavy-particle eigenstates i/>(p, 4>) now live on the 
potential energy surface U(p,<f>), and are eigenstates of 
the effective Hamiltonian (in the EAA): 




-1 

2/M 


/ ld_ d_ 1 d 2 A 
dpPdp ^ p 2 d(j> 2 ) 


U{P, <t>) + 


ap 


i<j 


( 10 ) 


where a = a/6/3/(3 + /3). In order to extract scatter¬ 
ing lengths we require a representation of the wavefunc- 
tion that gives the appropriate asymptotic cluster states. 
We write: if)(p, (j>) = t^=o P~ 1/2 fn(p)Xn(p; </>), where the 
channel functions Xn satisfy the hxed-p equation: 

^^ d g^ L + u (p^)X n = Un (p)X n ( 11 ) 


Due to identical particle symmetry, one only needs to 
consider the restricted range <j> £ [0, 7 t/6 ) with the bound¬ 
ary condition EDI: 


r 1 d X 
llm 

£^0 X d(j) 


pA(l + /J) 

( P ) 

7T / 6 — £ V2/3 

\Z + PJ 


1/6 


( 12 ) 


Finally, the four-particle problem is reduced to a set of 
coupled equations in p only: 

=l^ + QM + 2P M §-)f ip) (i3) 

+ u cff (p)/(p) = E EAA f(p) 


Here, U c ff is a diagonal matrix with elements U n (p) — 
I/8/I4 p , Pmn{P ) = {Xm\Xn)(f, and QmniP ) 

(Xm | Xn)tj>- The Born-Oppenheimer potentials obtained 
from Eqs. ([5| and @ are independent of A, and only 
need to be calculated once for a given /3. The A- 
dependence only appears through the boundary condi¬ 
tions 0 and ©• 

In Table [IJ we sumarize the boundary conditions on 
the channel functions Xn(/ , ;</)- Note that starting with 
positive parity noninteracting bosons, one achieves the 
infinitely repulsive “fermionized” limit as A —► —00 (since 
5 hl < 0). Equation (12) then gives a boundary condi¬ 
tion identical to that of noninteracting fermions, but of 
negative parity. Due to Pauli exclusion, identical (spin- 
polarized) fermions are insensitive to the zero-range in¬ 
teraction, and the boundary conditions for fermions are 
unaffected by A. 


III. NUMERICAL SOLUTIONS 


TABLE I. The boundary conditions on the channel functions 
Xn{p', <t>) for the non-interacting (NI), infinitely repulsive (IR) 
cases are summarized here. 



(+) parity 

(-) parity 


1SJ„- 0 


NI bosons 

[sL„«-° 

x\ l j >=n /6 = 0 


lC7 =0 

xl ^> = o = 0 

IR bosons 

x\(j, = 7 r/6 = 0 

xl<£— tt/6 = 0 


Nl q(, = 0 = 0 


NI fermions 

Xl^n/6 = 0 

xl 4> —7t/6 = 0 


for the Li-Cs mass ratio /3 -1 / 2 = 4.7, and A —> — 00 . 
Panel (a) shows the hyperradial potential curves U n (p) 
while panel (b) shows the effective H-H interaction u(x). 
Trimer (HHL) and tetramer (HHHL) bound states are 
indicated by dashed red lines in panel (a) and (b), re¬ 
spectively. Panels (b-c) and (d-e) are identical to (a-b), 
but for A = —1 and A = 0, respectively. Note that the 
potential curves U n (p) depend on A through the bound¬ 
ary condition ( |T2| . While u{x) is independent of A, the 
trimer binding energies depend on A through Eq. ([ 7 ]). 

Note that in the limit p —> 00 , the lowest potentials 
U n (p) approach the trimer energies, asymptotically rep¬ 
resenting atom-trimer channels (H+HHL). In the low¬ 
est atom-trimer channel we clearly see the appearance 
of a short-range potential energy barrier as A —> — 00 . 
This barrier is a direct result of the boundary condi¬ 
tion (12), and represents the repulsive effect of fermion- 
ization. We also see a set of potentials asymptotically 
approaching the dimer energy U n /Bi = — 1. These rep¬ 
resent the three-body (H+H+HL) channels. Three-body 
recombination at threshold is controlled by the lowest 
such potential, which in the large p limit behaves as 
U n (p) —► K^ lin /(2/44p 2 ). A full hyperspherical calcualtion 
would also produce a set of potentials approaching the 
zero-energy threshold asymptotically representing four 
free atoms (H+H+H+L). Our calculation cannot capture 
these potentials because the light particle is required to 
be bound. 

Numerical solutions to Eqs. 0, ( |TT| ) and ( p~3|) are 
found by expressing the desired wave functions as a sum 
over b-splines and solving the resulting generalized eigen¬ 
value problem. We extract q,ad and a at by matching 
the solutions ip(x) and fo(p) to a phase-shifted cosine in 
the asymptotic region, and extrapolating the zero-energy 
limit of k tan 6 1 /a (see ; 3D] for details.) 


Born-Oppenheimer and hyperradial potential curves 
provide a view of the few-body energy landscape that aids 
in the interpretation subsequent calculations. In Fig. B a - 
b), we show potential curves and bound state energies 


A. Scattering lengths bad and 4 at 

In Fig. [2j we show contour plots of the atom-dimer 
scattering length oad and the atom-trimer scattering 
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FIG. 1. (color online) Panels (a), (c), and (d) show hyper- 
radial potential energy curves U u (p) determined by solving 
Eq. (11i for A —> —oo, A = —1 and A = 0, respectively. Pan¬ 
els (b), (c) and (d) show the corresponding potentials u(x) of 
Eq. q4p. The dashed red lines indicate the energies of trimer 
and tetramer bound states. 


length oat on the plane formed by tan -1 (—A) and /3 -1 / 2 . 
We consider only the lowest atom-trimer channel for 
the elastic collision H + HHL —► H + HHL when cal¬ 
culating oat- All higher scattering channels are en¬ 
ergetically closed for all p , and their effect is negligi¬ 
ble in comparison to the EAA. The Born-Oppenheimer 
method requires the light particle to be bound such that 
the lowest atom-dimer scattering channel asymptotically 
contains an HL dimer, describing the elastic collision: 
H+HL —> H+HL. Our calculation fails below the dot¬ 
ted line in Figs. [2j a) and (b), which is described by the 
formula A = + f3). Below this line, the lowest 

atonr-dimer scattering threshold contains an HH dimer 
and a free L atom, and the structure of the three-body 
phase diagram is fundamentally different [29]. 

The dashed lines in Figs. 2](a) and (b) indicate oad —> 
0 and a at —► 0 , respectively, while the solid lines indi¬ 
cate |(jad| —t oo and IoatI —> oo. All points along the 
dashed lines are characterized by reflective elastic colli¬ 
sions in the threshold limit. Points along the solid lines 
are characterized by reflectionless collisions and the ap¬ 
pearance of a new bound state. The number of trimer or 
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FIG. 2. (color online) Panels (a) and (b) show contour 
plots of oad and hat, r espective ly, on the plane formed by 
tan -1 (<?hh/|<?hl|) and \Jmn/m l. The dashed lines mark the 
locus of points where the corresponding scattering length van¬ 
ishes, while the solid lines indicate where it diverges. Along 
the dotted line, the (HH) binding energy is equal to the (HL) 
binding. 
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FIG. 3. (color online) The A-dependence of various energies 
involving HL, HHL, and HHHL clusters for /3~ 1/2 = 4.7, cor¬ 
responding to the Li-Cs mass ratio. All energies shown are 
calculated in the EAA. 


tetramer bound states is indicated by n in each figure. 
The present calculation smoothly connects the poles and 
zeroes found along A = 0 and |A| —>■ oo in [50] . 

We see from Fig. [2] that the atom-dimer scattering 
length is typically much larger than the HL scattering 
length: oad a, while the same is not true of a at- We 
trace the reason for this back to the behavior of the po- 
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tential curves in Fig. |T] First, the trirner (ground-state) 
binding is deeper than the HL binding, and therefore the 
size of the trirner is correspondingly smaller because it 
is confined to the deeper portion of the potential u(x). 
This means that the effective atom-trimer interaction is 
of shorter range than the atom-dimer interaction, and 
typically a ad is also large in comparison to a at- 


B. Thresholds energies for few-body processes 

In Fig. |3j we show the three-body ground state energy 
£3 (blue dotted line) and the four-body ground-state en¬ 
ergy £4 (red dashed line), both calculated in the EAA 
as a function of tan -1 (—A)/ 7 r for the Li-Cs mass ratio 
/3 -1 / 2 = 4.7. As one might expect, the energies increase 
monotonically with increasing ^hh due to the increas¬ 
ing H-H repulsion. We also show the bound state en¬ 
ergy 2£ 2 for two dimers (HL+HL, black line), £ 2 + E 3 
for a dimer plus trirner (HL+HHL, black dot-dashed 
line), E 2 + £4 for a dimer plus tetramer (HL+HHHL, 
purple dot-dot-dashed line), and 2E 3 for two trimers 
(HHL+HHL, orange long-dash line). Note that for this 
particular mass ratio, when tan -1 (—A) < 0.147T, one 
finds E 3 < 2E 2 , while for tan -1 (—A) > 0.147T, E 3 > 2E 2 . 
The critical value tan -1 (—A) = 0.147T (marked by a 
black circle) indicates the transition point such that in 
a gas of dimers, trirner production through the process 
HL + HL —► HHL + L is energetically allowed for weaker 
H-H repulsion. Placing the critical value on the plane 
{/3 -1 / 2 , tan -1 (—A)/ 7 r} = {4.7,0.14} (see Fig. [4j) and re¬ 
peating the calculation for other mass ratios allows the 
parameter space to be partitioned into regions where this 
process ((10 in Table [IT} is either allowed or disallowed. 
We are able to estimate the error incurred by the Born- 
Oppenheimer approximation by numerically calculating 
the line E 3 = 2E 2 both in the EAA and the UAA, drawn 
as the bottom and top black lines, respectively, bounding 
the shaded band in Fig. [4] We find that the critical mass 
ratio in the |A| —> 00 limit is /l -1 / 2 ss 7.05 in the UAA 
and /3 -1 / 2 ss 7.11 in the EAA, bracketing the previously 
quoted value /3 -1 / 2 = 7.0593 [251 . 

One can imagine other inelastic processes (listed in Ta¬ 
ble 0 involving three, four, five and even six particles 
that have thresholds determined purely by the HL, HHL 
and HHHL binding energies. For instance, the rearrange¬ 
ment process (Fj 6 j) becomes energetically allowed when 
£4 < £3 + E 2 . For the specific mass ratio /l -1 / 2 = 4.7, 
£4 = £3 + £ 2 at tan (—A) = 0 . 087 T as indicated by the 
black square on Fig. [ 3 ] For arbitrary values of /3 -1 / 2 > 1, 
process (F§ becomes energetically allowed below the 
dashed blue curve in Fig. |4j The three-body recombi¬ 
nation processes ( 10 .( 1 ®, and ( 10 ) are nearly always 
energetically allowed for /W 1 / 2 > 1 because there is al¬ 
ways an available dimer, trirner and tetramer state. The 
only exception is a small region near the Tonks-Girardeau 
limit (|A| — > 00 ) where the first trirner state appears ex¬ 
actly at /3 -1 / 2 = 1 [29] (The Born-Oppenheimer approx- 



V m H / m L 


FIG. 4. (color online) The two-dimensional parameter space is 
partitioned by lines corresponding to the degeneracy of scat¬ 
tering thresholds. Those processes from Table [IT] that are 
energetically allowed are listed in each region. 


imation gives /? 1 / 2 ss 1.08 EL- and the first tetramer 
state does not appear until ,0 -1 / 2 > 1.4 ^Oj. Therefore, 
there is a small region in the top left corner of Fig. [4] 
(bounded by the solid green line) where process (F§ is 
not allowed because no tetramer state exists. The colli¬ 
sion of two trimers may produce a tetramer and a dimer 
through process (10 only for values of to the left 
of the star in Fig. [3j and below the orange dot-dot-dash 
line in Fig. [0 For the range of mass ratios considered 
here, (10 is not allowed for fermionic H atoms. Inelastic 
processes such as these can lead to atom loss and thermal 
heating of the trapped gas. Conversely, atom loss rates 
can be measured as a signature of such processes. 

We have labeled each region in Fig. [4] by the set of re¬ 
actions listed in Table [0 that are energetically allowed in 
the forward direction. At sufficiently low temperatures, 
the energy dependence of the reaction rates for these pro¬ 
cesses is governed by the corresponding Wigner threshold 
law. For each of the three-body recombination processes 
(010, the threshold law gives |<Sfi| 2 oc E^ in , where 
ftmin ^ The hyperangular quantum number for the lowest 
three-body hyperspherical harmonic in the limit of large 
hyperradius R —> 00 [3S] • In general, K m i n is an irrational 
number determined purely by the masses of the collision 
partners. We list the corresponding values of K m i n in Ta¬ 
ble [0 The energy dependence of scattering probabilities 
of two-body processes (1010 is controlled not by the 
momentum in the entrance channel, but rather in the 
fragmentation channel [351EZ]: I'Sfil 2 ^ &/, and therefore 
approaches a constant at threshold if the fragmentation 
channel lies below the entrance channel. With the excep¬ 
tion of (F§, each of the processes in Table 0 becomes 
allowed in the reverse direction on the opposite side of 
the corresponding critical line in Fig. [4j 
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Inelastic process 

Threshold law: 

is fl r oc icr 

(PI) H + H + L—>HL + H 

Kmin = 7T/I0 “ 7r / 2 
tan = y/P/(2 + P) 

(P2) HL + H + H HHL + H 

^min — 7 t/(2|0|) 

tan (j>= v / (l + /3)/(3 + /3) 

(P3) HHL + H + H -)■ HHHL + H 

^min — (21^1) 

tan 4> = \/(2 + /?)/(4 + P) 

(P4) HL+HL+H -> HHHL +L 

Kmin = n/\(f) — 7t/2 

tan^ = 1/^/3 + 2/3 

(P5) HL+HL ->• HHL + L 

^min == 0 

(P6) HHL + HL —¥ HHHL + L 

^min 0 

(P7) HHL + HHL -» HHHL + HL 

K'min == 0 


TABLE II. Various inelastic processes that are energetically 
allowed in the regions shown in Fig. [l] are listed here. 

IV. CONCLUSION 

We have calculated atom-dimer and atom-trimer scat¬ 
tering lengths, as well as trimer and tetramer energies 


for the homogeneous ID HHL and HHHL systems as 
a function of the mass ratio /3 = mj^/mn and cou¬ 
pling ratio A = 5hh/shl- We expect our purely ID 
calculation to be relevant to current experiments with 
atomic mixtures in 2D optical lattices in the “BCS” 
limit, when aj_/a 3 D —> — oo such that the ID dimer is 
weakly bound and has large axial extent compared to 
a±. When the few-body system contains only one L 
impurity, or one ignores the L-L interaction, all eigen¬ 
states are parameterized by the ratios A = <7hh / 3hl and 
/3 = TOi/m#. We have determined the regions in this pa¬ 
rameter space where certain inelastic processes involving 
HL, HHL and HHHL clusters are energetically allowed, 
potentially leading to atom loss or heating. 

Finally, we note that the potential curves shown in 
Fig- 0 can be used to calculate the energy-dependent 
scattering cross section for the three-body recombination 
process H+H+HL—^H+HHL. Such calculations are be¬ 
yond the scope of this paper, but may be pursued in the 
future. 
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thanks the Institute for Nuclear Theory at the University 
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